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Abstract 

We show that the dispersion relation for gradient-drift and Farley- 
Buneman instabihties within the approximation of the two-fluid MHD 
should contain the terms which are traditionally supposed to be small. 
These terms are caused by taking into account divergence of particles 
velocity and electron density gradient along the magnetic field direction. 
It is shown that at heights below 115km the solution of the dispersion 
relation transforms into standard one, except the situations, when the 
electron density gradient is parallel to magnetic field or wave-vector. In 
these cases the traditionally neglected summands to the growth rate of 
the irregularities becomes significant. 

The additional terms depend on relative directions of electron density 
gradient, magnetic field and mean velocities. This leads to the difi'erent 
instability growth conditions at equatorial and high-latitude regions of the 
ionosphere. The obtained results do not contradict with the experimental 
data. 

1 Introduction 

In the recent review |Farely,2009| the equatorial instabihties at E-layer heights 
was analyzed, and the fohowing problem was formulated: "How different are the 
auroral and equatorial zones? Is the physics essentially the same, except for the 
fact that the auroral zone electric fields are often much stronger? Or is it impor- 
tant that in the auroral zone there are gradients parallel to the magnetic field 
in electron density, electron and ion temperatures and collision frequencies?" . 
In this paper we try to answer to this question. 

As stated in |Farely,2009| , the analysis of gradient-drift instabilities in the 
ionospheric plasma at E-layer heights in linear approximation is based on the 
theory suggested in [Fejer et al.,1975| . The result of the work [Fejer et al.,1975| 
is the correct derival of the dispersion relation for the plasma instabilities in 
the specific case when electron density gradient is oriented vertically, magnetic 
field - horizontally, electrons are magnetized and ions are unmagnetized. In 



the work [Fejer et al.,1975 authors also retrieved the well known solution of 
the dispersion relation. The attempt of solving the wider problem, when elec- 
tron density gradients and magnetic field have arbitrary directions, was made 
in [Fejer et al.,1984| and is thought to be valid both for equatorial and high- 
latitude ionosphere. Unfortunately in the both papers the solution of the dis- 
persion relation was obtained under a number of additional assumptions. Due 
to this fact these theories sometimes do not absolutely correctly describe the 
conditions of the instabilities generation, especially when the instabilities wave- 
vector or magnetic field are almost parallel to electron density gradient, as 
mentioned in |Farely,2009| . The possible cause of these problems is preliminary 
simplification of the dispersion relation and neglection of important terms. 

The goal of the paper is the analysis of the dispersion relation solution, gen- 
eralizing the formulas obtained in [Fejer et al.,1975 Fejer et al.,1984[ by taking 
into account some terms traditionally neglected. We base on the same two-fluid 
magnetohydrodynamics (MHD) approximation as used in [Fejer et al.,1975| and 
assume that this approximation is valid for the heights under investigation. We 
use approach suggested in |Gershman,1974[ |Gelberg,1986 for studying the polar 
ionosphere. Then we investigate all the terms in the obtained solution to take 
into account important ones. 



2 Solution of dispersion relation 
2.1 Traditonal approach 

The instability growth condition in presence of arbitrary direction of electron 
density gradient and electric and magnetic fields in linear approximation is 
retrieved, for example, in [Fejer et al.,1984| . The dispersion relation obtained 
there has a form: 
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where = Veo — V^q - relative electron/ion drift velocity; Cg — y "^"m 
ionacoustic speed, a - recombination coefficient; 
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a coefficient defining the aspect sensitivity of the instabilities, b = ^. 
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To emphasize the problem of the dispersion relation Q, let us analyze a 
special case, when electron density gradient is parallel to the wave-vector and 
perpendicular to the magnetic field. In this case we have 

k_L- (V7V X b) = (3) 

and in the standard dispersion relation ([I]) the presence of gradients can be 
neglected, and gradient-drift generation mechanism stop working. 

The condition ^ is almost valid for non-disturbed equatorial ionosphere 
(where electron density gradient and wave vector are almost vertical), and for 
non-disturbed high-latitude ionosphere (where electron density gradient and 
magnetic field are almost vertical), so this problematic situation should arise 
frequently. 

But sometimes this contradicts with the experimental data, for example at 
Earth magnetic equator, with vertical gradient of electron density and vertical 
wave-vector of sounding wave (and vertical wave-vector of irregularities, corre- 
spondingly) a powerful scattering is observed, as stated in |Farely,2d09l pp. 1515- 
1516]. This allows us to suggest the presence of the gradient-drift instabilities 
even in this special case. 

It is obvious that in the process of the dispersion relation ([T]) output a number 
of approximations was used, and necessary effect may be lost or neglected due 
to its small value. 

As a basis for our analysis we will use the approximation of two-fluid MHD, 
following to the traditional approach developed in [Fejer et al.,1975 Fejer et al.,1984 



But in opposite to this approach, following to the works Gershman,1974 Gelberg,1986| 



we will not simplify the dispersion relation before obtaining the solution. 
2.2 Full solution 

Usually, when investigating ionospheric instabilities, for example, in |Farely,2009| 
it is thought that gradient-drift instabilities should be investigated when ions are 
unmagnetized, and electrons are magnetized, and velocity divergence effect can 
be neglected (Rogister and D'Angclo,19701 . But it is well known that in fact the 
dispersion relation in two- fluid MHD approximation has more general form, sim- 
ilar to the one obtained in [Fejer et al.,1975| , but with taking into account the 
case of arbitrary magnetized particles of both types (i.e., without preliminary 
simplifications of the dispersion relation) and average velocity divergence (see, 
for example, high-latitude case investigated in Gershman,1974[ [Gelberg, 1986| ) . 



It can be shown (see, for example, approach described in Gelberg, 1986[ 
sect. 5]) that in this case the solution of the dispersion equation will be more 
complex: 



VeO+W,o, ,. * 1 //Voe+*Voi,V ^2,2' 
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where = V^o — V^g - drift velocity of the electrons relative to ions. The 
solution was obtained from system of two-fluid MHD equations for both types 
of particles (ions and electrons) , according to the consistency condition for the 
system, in long wavelength approximation allowing to consider plasma density 
variations as quasi-neutral ones {6Ne — SNi). 

It is important to note that average velocities V^q i^i traditional approxima- 
tions of equal ionization and recombination Q — aN^Q = and smooth velocity 
profile (VqoV) = (see, for example, [Golant et al.,1980] ) can be defined 
as 

where CTq,, Dq, D^q are the the operators of conductivity, diffusion and ther- 
modiffusion correspondingly, discussed, for example, in [Golant et al.,1980] , and 
^ct'n = m "+rn ^an " IS a uormcd frequency of elastic collisions with neutrals. 

In case, when Q — aN^g ^ and/or (V„oV) ^ 0, the velocities are 
defined by solving the following well known zero-order equations numerically or 
analytically: 



Z^eNaoEo + Z^eNaoVao x ^0+ 

+m„iV„o(V„o-V)VaO + ViTaNao) + iV„oma(V„o - VnV^i^ = ^ ' 

= Z,Nm + NeO 

It is important to note, that when investigating F-B and G-D instabilities 
in two-fluid MHD approximation usually the following condition is expected to 
be valid: 

\Q-aNl,\«aN^^„Q (7) 
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If that not the case, we should take into account densities of the differ- 
ent kinds of particles (electrons and ions) in the ionization/recombination term, 
which usually leads to the ionization waves, discussed, for example, in [Akhiezer et al,1974] . 



3 Discussion 



3.1 Full solution 



By taking into account only significant terms at the E-layer altitudes, we obtain 
from equation Q the following simplified solution for FB and GD instability: 




~^ {2aNo) + ^J^f^ ) ^ ( ^ ) . (b x k) 
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+.7^^ff^-b)(b-k)| (8) 



'{l + 'f)^mA'J^k^ \\ No 
The first five summands are well known: 

first summand defines Farley-Buneman instability frequency in form, ob- 
tained, for example, in [Fejer et al.,1975| ; 

second one defines the growth rate of Farley-Buncman instability derived in 
[Farley, 1963[IBuneman T963] : 

third one defines the affect of recombination processes to the growth rate, 
obtained, for example, in [Fejer et al.,1975| ; 

forth one defines the effect of electron density gradients due to their arbi- 
trary orientation to the magnetic field and wave-vector, according to the paper 
[Fejer et al.,1984[ . In a simple case of gradients perpendicular to the magnetic 
field it transforms to the widely used (as stated in [Farely,2009| ) summand ob- 
tained in [Fejer et al.,1975[ . 

fifth one defines velocity convergence effect. Usually the term is neglected. 
It is well known that the term actually can be neglected in case of gradient-free 
case IRogister and D'Angelo,1970|, but in general case it should be taken into 



account, see, for example, Gershman, 1974[ , [Gelberg, 1986| eq.5.15] 



In the equatorial ionosphere the last term can be neglected due to the mag- 
netic field and electron density gradients are almost orthogonal and the term is 
small. But in general case, valid for both polar and equatorial ionosphere, one 
should take it into account. 

It is important to note that the first four summands in the solution are 
well known, and analyzed, for example, in [Fejer et al.,1984[ , but the last two 
summands affects only on growth rate and becomes significant in the regions 
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where electron density gradient is parallel to the magnetic field or wave-vector, 
and traditional gradient term, obtained in [Fejer et al.,1975 Fejer et al.,1984| 



(forth summand in the expression ^ ) becomes zero. 

It is also important to note that the second summand, which defines F-B 
mechanism of instability growth, decrease exponentially with altitude due to 
decreasing the ratio jt^^tttt with increasing the height. From the other side, 

the fifth term does not decrease so quickly so it depends mostly from the relative 
gradient of electron density . This leads to the importance of the fifth term 
at higher altitudes. 

As one can see, the forth term, which defines traditional G-D decrement, 

also decrease exponentially with altitude due to decreasing the ratio (^-ff-^ ■ At 
the same time, the sixth term does not have such a fast changes with altitude, 
due to — 't'", , relation remains almost constant with altitude. So the sixth 
term becomes also important at higher altitudes (above 100km). 

Let us qualitatively analyze the obtained solution ^ by supposing that 
electron density gradient is almost vertical. This situation corresponds to the 
non-disturbed E-layer, produced mostly by ionization/recombination processes. 
By doing this we neglect the presence of turbulence and wave-like processes, 
which should be taking into account during more detailed analysis. 

3.2 Equatorial case 

At the equator we can neglect the sixth term in ([8|, and the solution becomes 
a bit simpler: 



1 * 1 + * z.*'^ I V 1 + * 



Let us analyze the situation, when magnetic field direction b is horizontal, 
wavevector k is nearly vertical, vertical average velocity Voe_L = Voi_L = 20m/ 
horizontal electron/ion drift velocity at 100-110 km Vd = Vbe — Voi — 200m/s 
(the geometry and physics of this is discussed, for example, in |Farely,2d09] ), 
vertical gradient is supposed to be 10^''m~^, k — 0.1m~^and almost verti- 
cal. It is important to note, that in equatorial case the electron/ion drift velocity 
Vd is strictly horizontal (see, for example, |Farely,2009| ). From the other side, 
the average velocity has a vertical component, and the last term in Q is always 
not equal to zero. 

At Fig.l we show a height dependence of the second, forth and fifth terms in 
([9| on the altitude for the described conditions. Geometry is shown at Fig. 1(11). 
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Figures IB-D corresponds to different zenith angle in east-west plane (flow angle 
/3, correspondingly 1.0,10.0,30.0 degrees, for zero aspect angle), figures E-G 
corresponds to the different zenith angle in north-south plane (aspect angle a, 
correspondingly 0.1,1.0,3.0 degrees, for zero flow angle). 

Figure lA corresponds to both angles equal to zero. As one can see, at all 
heights the velocity divergence term '^""i^X'"'^ '^^ (dot-dashed line) in this 

case is much larger than standard F-B increment term j^-jhi ^"l+'tti^"' 
(dotted line) and due to this must be always taken into account. At heights 
above 115-120 km it becomes larger than decrement term 'i^^s^^ (solid 

line) and this fact can cause an additional growth of the instabilities. This term 
can be important, for example, for analysis of the echo at heights above 115km 
(see, for example, [Farely, 20091 Fig-3])- 

As one can also see, at flow angles higher than 10 degrees (Figure IC-D) 

standard GD term ^^''^'^2 ("ff^) { (^^^) ' ^ ^)} (dashed line) becomes im- 
portant for generation of the instabilities, and it is well known fact (see, for 
example, |Farely,2009| ). Also, as one can see, the effect of the new term will be 
stronger at larger wavelengths and weaker at smaller wavelengths. 

It can be easily shown, that for lower part of ionospheric E-region, when 
electron density gradient and velocities are co-directed (this corresponds to the 
day conditions of the equatorial ionosphere) , the last term increases the growth 
rate. In this case the possibility of observation of the irregularities becomes 
higher. When the electron density gradient and velocities are anti-directed (this 
corresponds to the night conditions of the lower part of E-layer of the equatorial 
ionosphere), the last term decreases the growth rate. In this case the possibility 
of observation of the irregularities becomes lower. 

This effect - increasing of irregularities level at day and decreasing at night 
is observed at the 100-110 km. altitudes more-less regularly (see, for example, 
in |Farley,1985] |Abdu et al.,2002| [Denardini et al.,2005| ). The similar hourly 
dependence of the irregularities is observed even at higher altitudes at equator 
(see, for example, in |Chau and Kudeki,2006| ). It is important to note that 
daytime and nighttime power dynamics at 100-1 10km is explained for non- 
perpendicular wavevector and mean electron velocity by equatorial electrojet 
dynamics (see, for example, in |Farely,2009| ), similar dynamics of scattered 
power at higher altitudes (above 130km) is still unexplained (see, for exam- 
ple, in Chau and Kudeki,2006| ). So the new term have similar dynamics and 
does not contradict with these known experimental results. 



3.3 Quasi-polar case 

Let us analyze the following case: strong geomagnetic disturbance produces 
almost horizontal Voe = Vd = 400m/s,Voi = 0, electron density gradient 
and magnetic field are nearly parallel (within 10 degrees) and almost vertical, 
wavevector nearly perpendicular to the magnetic field, — 10~''m~-'^and 

vertical, wavevector k = O.lm^^ and almost horizontal. In this case we can 
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see that the last summand in the solution ([s]) is significant and the solution (|8| 
becomes: 
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ujr + iT = T^•k + ^— I ( ,^-k) -Cffc2 




(b X k) 



(10) 



As one can see, the last summand in expression ( 10 ) can produce instabilities 
at aspect angles larger than traditional theory predicts. The "large aspect 
angles" effect is known in high latitude ionosphere and has different explanations 
(see, for example, in [Hamza and St. Maurice, 1995| and references there). As 
also one can see, the new term do not depend on wavenumber, and standard 
F-B term increases with wavenumber growth. Due to this the presence of the 
last term will be more significant at larger wavelengths than at shorter ones. 

At Fig. 2 the numerical results are shown. Electron density gradient is sup- 
posed to be vertical, magnetic field does not have east-west component and 
is about 10 degrees off the zenith. Maximal drift velocity is 400m/s and ori- 
ented in north-south direction, which corresponds to average geomagnetic dis- 
turbance. Geometry is shown at Fig. 2(H). Fig.2B-D shows dependence in the 
angle in plane perpendicular to the magnetic field (flow angle /3, 1.0,3.0,10.0 
degrees correspondingly, aspect angle 1.0 degree), Fig.2E-G shows dependence 
in plane of magnetic field line (aspect angle a, 1.0,3.0,5.0 degrees correspond- 
ingly, zero flow angle). Fig 2 A corresponds to the zero flow and aspect angles. 
As one can see, the new term (dot-dashed line) is not important at zero as- 
pect angles, but becomes important at aspect angles 1-3 degrees (Fig.2E-F), 
and at these aspect angles can be responsible for generation of instabilities. 



when standard F-B increment term ^"'14-'^^°' (dotted line) is 

smaller than decrement term j^-^C^k^ (solid line). Due to specific ori- 
entation of the electron density gradient and magnetic field, the standard G-D 
term l"^) {(^1^) ' ^ (hashed line), is almost unimportant at 

these small flow angles. 

As one can see, the new term can cause generation of the instabilities at 
larger aspect angles (Figs.2E-F). 



4 Conclusion 

Our analysis has shown that growth rate for gradient-drift instabilities has 
significant terms that depends on relative direction of electron density gradi- 
ents, mean velocity and magnetic field directions pi. This leads to the dif- 
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ferent conditions for the growth of the F-B and G-D instabihties in the equa- 
torial and auroral zones of the ionosphere. These terms are usually neglected 
IFejer et al., 19751 |Fejer et al.,1984|rFarely,2009| , but for building a general the- 
ory they must be taken into account. 

At heights below 115km the full solution (|8| transforms into standard one 
[Fejer et al.,1984| , except the situations, when the electron density gradient is 
parallel to magnetic field or wave-vector. In these cases the new terms in the 
growth rate of the instabilities becomes significant and predicts the presence of 
gradient-drift instabilities even in these cases, as shown by expression ([8| , in op- 
posite to traditional theory obtained in [Fejer et al.,1975| [Fejer et al.,1984| . At 
low-latitude region these terms may cause additional growth of the instabilities 
at higher altitudes (above 115-120km). At high-latitude region these terms may 
produce changes of irregularities aspect sensitivity (large aspect angles effect). 
Both effects do not contradict with experimental data. 
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Figure 1: The absolute value of each term in the solution ([9| in the equatorial 
ionosphere for different aspect (a) and flow {(3) angles. Figures (B-D) corre- 
sponds to different flow angles, correspondingly 1.0,10.0,30.0 degrees, for zero as- 
pect angle. Figures (E-G) corresponds to the different aspect angles, correspond- 
ingly 0.1,1.0,3.0 degrees, for zero flow an^le. Figure (A) corresponds to both an- 
gles are equal to zero. Terms: Y^T^C^fc^ - solid hue; j^-jitf ( ^"l+^j^"' 

- dotted line; (^) { j • (b x k)} - dashed line; M^.S^ _ 

dot-dashed line. At Fig. (H) the geometry and explanation of a and /? is shown. 
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Figure 2: The absolute value of each term in the solution (10) in the high- 
latitude ionosphere for different aspect (a) and flow (/3) angles. Figures (B-D) 
corresponds to different flow angles , correspondingly 1.0,3.0,10.0 degrees, for 1.0 
degree aspect angle. Figures (E-G) corresponds to the different aspect angles, 
correspondingly 1.0,3.0,5.0 degrees, for zero flow angle. Figure (A) corresponds 
to both 1.0 degree flow angle and zero aL^ect angle. Terms: ^r^iC^k'^ - solid 

b ) (b-k)) - dot-dashed line. At Fig. (H) 



line; 
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the geometry and explanation of a and /? is shown. 



